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An  experimental  investigation  was  performed  to  determine, 
in  situ,  the  change  in  mean  particle  size  across  the  exhaust 
nozzle  of  a  small  metallized  solid-propellant  rocket  motor. 
Light  scattering  profiles  were  recorded  at  both  the  exhaust 
and  the  entrance  of  the  nozzle.  The  experimental  method 
used  provided  excellent  results  within  the  exhaust.  However, 
combustion  light  at  the  wavelength  of  the  transmitted  light 
hampered  light  scattering  measurements  within  the  motor. 
Particle  size  measurements  were  consistent  with  the  sizes 
found  in  the  collected  exhuast  products. 
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I.  INTRODUCTION 

The  importance  of  aluminum  oxide  (AJ^^)  particle  size 
produced  in  solid-propellant  rocket  motors  has  been  recog¬ 
nized  since  aluminum  was  first  introduced  as  a  fuel  additive. 
Aluminum  added  to  a  propellant  will  increase  the  thrust 
produced  and  aid  in  suppressing  high-frequency  combustion 
instability.  Although  the  performance  improves  with  aluminum 
addition,  two  phase  flow  losses  adversely  affect  specific 
impulse  efficiency.  Specific  impulse  efficiency  is  the  product 
of  several  component  efficiencies,  i.e.,  chemical-kinetic, 
combustion,  two-phase  flow,  two-dimensional  flow,  submergence 
and  throat  erosion.  In  aluminized  propellants,  the  two-phase 
flow  loss  is  one  of  the  more  significant  losses.  Two-phase 
flow  losses  result  from  particle-gas  thermal  and  velocity  lags. 
The  latter  strongly  depends  on  particle  size.  Parametric 
studies  involving  small  motors,  with  throat  diameters  less 
than  one  inch,  at  AFRPL  [Ref.  1]  and  CSD  [Ref.  2]  have 
revealed  that  a  1  micron  change  in  particle  diameter  results 
in  a  change  of  1%  in  Isp.  The  effective  use  of  particle 
damping  for  suppression  of  combustion  pressure  oscillations 
is  also  strongly  dependent  upon  the  particle  size  distribution. 

One  of  several  analytical  computer  programs  used  to  pre¬ 
dict  performance  of  a  solid-propellant  rocket  motor  is  the 
SPP  (Solid  Performance  Program)  [Refs.  3,4].  As  the  above 


discussion  indicates,  it  is  of  critical  importance  to  accurately 
predict  particle  size.  This  is  especially  true  for  particle 
damping  with  minimum  smoke  propellants  which  typically  con¬ 
tain  less  than  2%  metal  by  weight.  The  1975  version  of  the 
SPP  [Ref.  3]  contains  a  particle  sizing  empirical  model 
developed  by  Cohen.  This  model  was  based  on  correlating 
collected  exhaust  samples.  After  a  critical  review,  this  model 
was  considered  inadequate.  A  new  model  was  therefore  required 
for  a  new  SPP  version,  titled  Improved  SPP  [Ref.  4] .  Inorder 
to  meet  SPP  accuracy  goals,  this  model  had  to  predict  particle 
size  within  ±10%.  Unfortunately,  there  are  no  adequate 
theoretical  models  relating  particle  size  to  propellant  and 
motor  parameters.  This  necessitates  the  use  of  empirical 
methods  for  correlating  particle  size. 

There  are  two  methods  that  can  be  employed  to  develop  an 
empirical  relation.  One  is  to  correlate  particle  size  data  to 
propellant  and  motor  variables,  as  was  done  by  Cohen.  The 
other  method  is  to  calculate  particle  size  based  on  a  critical 
Weber  number  for  a  maximum  stable  drop  size  [Ref.  5].  After 
considering  the  relative  merits  of  each  method,  Hermsen  [Ref. 

5]  devised  an  empirical  model  based  on  collected  particle  size 
data  for  the  improved  SPP.  This  model  was  correlated  to  a  much 
broader  data  base  than  was  previously  available  to  Cohen. 

The  particle  diameter  correlated  was  the  mass-weighted 


diameter,  given  by: 
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Smith  [Ref.  6]  has  shown  that  was  the  most  appropriate 
average  to  use  in  predicting  two-phase  flow  losses. 

Hermsen's  empirical  relation  [Ref.  5]  is  presented 
below: 


D43  =  3.6404  D;2932  (1  .e-.00081635PT)  _ 


where : 

D4;j  =  mass-weighted  average  diameter,  pm 
Dt  =  nozzle  throat  diameter,  in. 

£  -  Ai.2°3  concentration  in  chamber,  g-mol/100  g 

P  =  chamber  pressure,  psia 

t  =  average  residence  time,  msec 

The  predicted  D43  is  now  applied  to  the  expression  below  in 
order  to  calculate  the  two  phase  flow  loss  [Ref.  7] . 
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where : 


ETATP  =  impulse  loss  effect,  % 

£  =  mole  fraction  of  condensed  species,  moles 

of  condensed  species  per  100  grams  of 
mixture 

=  D.,  in  microns 

p  4  j 


P  =  average  motor  chamber  pressure,  psia 

e  =  nozzle  expansion  ratio  at  ignition 
conditions 

Dfc  =  nozzle  throat  diameter  at  ignitiion, 
inches 

e  =  Naperian  base  2.71828 

C(  )  =  coefficients  that  are  a  function  of  D 

and  Dp 

The  two-phase  flow  loss  is  applied  with  the  other  aforementioned 
efficiencies  to  a  theoretical  Isp.  The  Isp  loss  from  boundary 
layer  effects  is  then  subtracted  from  the  corrected  theoreti¬ 
cal  Isp  to  obtain  delivered  Isp. 

Collecting  exhaust  products  is  feasible  only  for  small 
motors  and  the  techniques  used  often  result  in  variation  in 
the  measured  size.  Validation  of  the  SPP  model  must  be  based 
on  data  collected  within  the  motor  environment.  Light 
scattering  measurements  are  well-suited  for  this  application. 
These  measurements  are  more  easily  made  if  the  particles  are 
large  compared  with  the  wave  length  of  the  transmitted  light 
and  the  transmittance  is  greater  than  90%.  The  former  allows 
the  accurate  application  of  Fraunhofer  diffraction,  whereas  the 
latter  specification  is  necessary  in  order  to  satisfy  single 
scattering  requirements . 

The  Naval  Postgraduate  School  has  conducted  a  series  of 
investigations  to  determine  the  applicability  of  light  scatter¬ 
ing  measurements  to  the  solid  rocket  motor  exhaust  products 


[Refs.  8,9,10,11].  This  investigation  was  directed  at  modi¬ 
fying  the  apparatus  and  techniques  developed  in  the  previous 
efforts  in  order  to  more  accurately  determine  particle  size. 
Validation  of  the  apparatus  was  also  a  goal.  The  SPP  particle 


II.  THEORETICAL  BACKGROUND 


The  most  general  theory  of  light  scattering  was  developed 
by  Mie  and  is  presented  by  Van  de  Hulst  [Ref.  12].  Mie's 
theory,  although  not  restrictive  in  its  assumptions  in  regard 
to  particle  size  and  refractive  index,  is  however  very  com¬ 
plex  to  utilize  for  data  reduction.  Dobbins,  et  al.  [Ref.  13] 
and  Powell,  et  al.  [Ref.  14]  derived  similar  expressions  for 
intensity  ratios  at  two  forward  scattering  angles  based  upon 
an  Upper  Limit  Distribution  Function  (ULDF) .  and  Fraunhofer 
diffraction.  The  ULDF  was  initially  proposed  by  Mugele  and 
Evans  in  Reference  15.  Although  the  expressions  developed  by 
Dobbins  and  Powell  are  extremely  useful  in  relating  intensity 
ratios  to  D^ •  the  volume-surface  mean  diameter,  these  models 
are  still  somewhat  cumbersome  to  use  for  data  reduction 
purposes . 

Particle  sizes  found  in  solid-propellant  rocket  motors 
are  often  bi-modal.  Particles  in  the  larger  mode  are  the  most 
important  for  two-phase  flow  loss  calculations  and  are  gener¬ 
ally  much  greater  in  diameter  than  the  wavelength  of  inci¬ 
dent  light.  Scattering  by  these  large  particles  is  described 
by  Fraunhofer  diffraction.  Hodkinson  [Ref.  16]  and  others 
[Refs.  13,14]  have  demonstrated  that  diffraction  theory 
(within  the  forward  lobe)  can  be  used  to  approximate  the  Mie 
theory.  Buchele  in  Reference  17  summarized  the  experimental 

15 


techniques  for  determining  the  particle  size  by  measurements 
of  dif fractively  scattered  light.  Buchele  presents  a  function 
which  closely  approximates  the  curves  given  in  References  13 
and  14. 


1(6)  =  EXP [- ( .57  a  8)2] 


where : 


1(0)  =  ratio  of  intensity  of  scattered  light  at 

some  angle  (6)  to  the  intensity  of 
scattered  light  at  theta  equal  to  zero 
degrees 

a  =  t\D/X  is  the  particle  size  parameter  for 
diameter  D,2  and  wavelength  of  light 
lambda  ( X ) 7 


The  above  expression  was  used  in  the  current  study  for  data 
reduction  to  obtain  D32.  Unfortunately,  D43  is  the  mean 
diameter  used  in  the  SPP.  Therefore,  D^2  must  be  related  to 
D43  either  by  analyzing  exhuast  samples  or  by  approximate 
analytical  means. 


A  final  important  point  is  that  the  above  expression  is 
valid  for  a  <_  3  (the  center  lobe)  .  For  the  expected  ranges 
of  theta  and  D,2,  this  proves  to  be  a  very  mild  restriction. 


III.  EXPERIMENTAL  APPARATUS 


A.  ROCKET  MOTOR 

The  rocket  motor  used  in  this  investigation  was  the  same 
as  the  long  motor  described  in  Reference  11.  The  motor  com¬ 
ponents  are  presented  in  a  photograph  in  Figure  3.1  and  the 
installed  motor  is  shown  in  Figure  3.2.  A  schematic  is  also 
presented  in  Figure  3.3.  The  propellant  grain  dimensions  re¬ 
mained  the  same  as  in  the  previous  effort  [Ref.  11] .  The 
grain  was  two  inches  in  diameter,  one  inch  in  length  and  had 
a  web  thickness  of  .725  inches.  The  grain  was  cylindrically 
perforated  and  allowed  to  burn  also  on  the  aft  end.  The 
igniter  was  composed  of  BKNO^ •  Ignitiion  of  the  igniter  was 
achieved  by  heating  a  resistance  wire  with  a  12  VDC  source. 

Some  of  the  exhaust  products  were  collected  in  an  eight 
inch  diameter  stainless  steel  exhaust  tube.  The  products 
were  collected  from  both  the  side  walls  of  the  tube  and  the 
end  cap.  Another  particle  sample  was  collected  from  the  walls 
of  the  combustor,  near  the  nozzle  entrance.  Both  samples  were 
then  cleaned  and  mounted  for  observation  using  a  scanning 
electron  microscope. 

B.  LIGHT  SCATTERING  APPARATUS 

The  apparatus  shown  in  Figures  3.4  and  3.5  was  basically 
the  same  as  used  in  Reference  11  with  some  exceptions.  The 


light  source  used  for  the  nozzle  exhaust  was  a  8-miliwatt, 
helium-neon  laser.  The  nozzle  entrance  also  used  a  8-milliwatt 
helium-neon  laser.  The  lasers  were  mounted  on  two  parallel 
optics  benches,  such  that  one  beam  could  pass  through  the  ex¬ 
haust  and  the  other  beam  could  pass  through  the  motor  cavity 
adjacent  to  the  nozzle  entrance.  A  beam  expander/collimator 
was  used  with  the  helium-neon  laser  in  the  exhaust.  Each  beam 
was  intercepted  by  a  stop  placed  directly  in  front  of  the 
narrow  pass  filter.  The  narrow  pass  filter  was  used  to  eliminate 
extraneous  light.  The  stop  and  narrow  pass  filter  were  placed 
30.5  centimeters  from  the  motor  centerline  for  the  exhaust 
beam  and  11.5  centimeters  from  the  motor  centerline  for  the 
motor  cavity  beam.  The  exhaust  beam  narrow  pass  filter  had 
a  diameter  of  5.08  centimeters  and  the  motor  cavity  narrow 
pass  filter  had  a  diameter  of  2.3  centimeters.  Both  beams 
passed  through  the  centers  of  these  filters.  In  the  present 
arrangement  of  the  apparatus  this  resulted  in  a  maximum  angle 
for  collecting  scattered  light  of  4.8  degrees  in  the  exhaust 
and  6.0  degrees  in  the  motor  cavity.  The  placement  of  the 
optics  closer  to  the  motor  centerline  allowed  for  a  greater 
maximum  angle  to  be  recorded  than  in  the  original  design. 

Directly  following  the  narrow  pass  filters  were  condensing 
lenses.  Both  paths  used  a  lens  with  a  diameter  of  5.08 
centimeters.  These  lenses  had  a  focal  length  of  50  centimeters 
and  were  utilized  to  focus  scattered  light  onto  a  linear  array 
of  photodiodes.  The  array  consisted  of  1024  silicon  photodiodes 


with  .025  millimeter  spacing.  The  apparatus  in  Reference  11 
positioned  the  lenses  such  that  only  angles  from  0  to  .04 
radians  could  be  observed.  In  actuality,  however  (as  a  conse¬ 
quence  of  diffraction  around  the  beam  stop  and  beam  spreading) , 
the  range  of  useful  angles  where  only  scattered  light  occurred 
was  from  .02  radians  to  .04  radians.  In  the  current  arrange¬ 
ment  the  diode  array  was  placed  1.0  centimeters  below  the  main 
beam  and  the  collecting  tube  diameter  was  enlarged.  This 
allowed  a  range  of  angles  from  .02  to  .07  radians  to  be 
observed. 

C.  DATA  ACQUISITION  AND  DATA  REDUCTION 

The  data  acquisition  system  is  presented  in  detail  in 
Reference  10.  The  system  controller  was  an  HP  9836S  computer 
and  the  HP  6942A  multiprogrammer  provided  A/D  conversion  and 
storage.  The  data  acquisition  system  was  programmed  to  take 
data  at  a  specified  point  in  the  motor  firing  by  inputing  a 
trigger  pressure  and  a  time  delay  after  the  trigger  pressure 
was  reached.  The  number  of  diode  scans  for  the  exhaust  and 
motor  cavity  were  eight  and  four  respectively  (limited  by 
the  current  data  acquisition  system  memory  size) . 

Two  methods  were  used  in  References  10  and  11  for  data 
reduction,  one  of  which  was  altered  for  this  study.  The  modi¬ 
fied  method  pertained  to  the  iterative  graphics  technique. 

In  References  10  and  11  the  experimental  data  were  normalized 
and  compared  to  a  theoretical  profile  that  was  also  normalized 
to  some  assumed  scattered  light  intensity  at  zero  degrees 


for  a  specified  032*  A  "best-fit"  of  the  experimental  profile 
to  the  theoretical  one  yielded  The  current  investigation 

determined  that  this  method  resulted  in  ambiguity  in  deter¬ 
mining  as  a  consequence  of  assuming  a  centerline  intensity 
for  the  scattering  profile.  The  method  used  in  the  present 
study  was  one  in  which  the  experimental  data  were  not  altered. 

A  theoretical  curve  was  generated  for  any  specified  value  of 
D^2  based  upon  a  selected  initial  angle  (0^)  and  the  correspond¬ 
ing  intensity  at  that  angle  ( X ^ )  from  the  unaltered  experimental 
data.  The  approximate  intensity  profile  for  a  polydispersion 
given  by  Reference  17  is: 


I2  =  I1EXP-{  (02  —  ej)  (.57  Dtt/A)2} 


Specification  of  a  series  of  values  for  02  results  in 
values  for  I2  and,  therefore,  an  I  vs  0  theoretical  profile 
that  passes  through  the  experimental  data  at  1^,  0]_.  If 
this  curve  does  not  fit  the  experimental  data,  then  other  values 
of  D^2  can  he  chosen  and  a  new  theoretical  curve  generated. 

A  "best-fit"  method  is  subsequently  used  to  determine 
In  summary,  the  method  used  for  obtaining  the  mean  particle 
size  is  based  upon  fitting  an  approximate  (although  quite 
accurate  within  specific  angle  restrictions)  theoretical 
profile  for  a  polydispersion  to  the  unaltered  profile  obtained 
experimentally  over  a  small  range  of  scattering  angles.  The 
modified  data  reduction  program  is  presented  in  the  Appendix. 


I 


The  second  method  was  similar  to  the  first  and  also  used 


the  above  expression  solved  for 


-Ln(I2/I1)  U/.57tt)2/(02  -  62) 


The  procedure  used  was  to  first  select  a  starting  angle,  9^, 
along  with  a  02/0^  rat^-°*  A  value  of  I2  at  02  was  obtained 
from  the  experimental  data.  This  was  repeated  with  increasing 
values  of  @2  and  then  again  with  different  values  of  02/0^. 

The  results  were  plotted  graphically  as  particle  size  vs 
scattering  angle  for  each  angle  ratio.  A  nearly  constant 
value  of  D32  vs  0  indicated  a  good  correlation. 
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Schematic  of  the  Motor 


Figure  3.5.  Schematic  Diagram  of  Light  Scattering  Apparatus 
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LIGHT  SOURCE  SPECIFICATION 


Helium-Neon  Laser  (Exhaust) 
1.  Manufacturer 


2.  Model: 


3 .  Type : 

4.  Output  Power: 


5.  Beam  Diameter: 


6.  Beam  Divergence: 


Spec tra-Phy sics 


Ke-Ne  Class  IIIB 


8  mWatt 


.92  mm 


.87  nr ad. 


Helium-Neon  Laser  (Motor  Cavity) 


1.  Manufacturer: 


2 .  Model : 


3 .  Type : 

4 .  Output  Power : 

5.  Beam  Diameter: 

6.  Beam  Divergence: 


Uni-Phase 


1305P 


He-Ne  Class  IIIB 


5  mWatt 


.  81  mm 


1.00  mrad. 


IV.  RESULT  AND  DISCUSSION 


A.  INTRODUCTION 

The  purpose  of  this  study  was  to  improve  and  validate 
the  accuracy  of  the  apparatus  used  for  particle  sizing  in 
the  rocket  motor  environment.  Also,  the  experimental  data 
will  provide  a  basis  to  verify  the  particle  size  prediction 
model  of  the  SPP. 

The  apparatus  had  been  calibrated  by  Harris  [Ref.  10] 
and  Kertadidjaja  [Ref.  11]  using  glass  and  polystyrene  spheres 
in  a  suspension  of  water.  Kertadidjaja  [Ref.  11]  also 
attempted  to  demonstrate  the  applicability  of  the  apparatus 
to  the  combustion  environment.  This  investigation  continued 
the  earlier  studies  with  several  refinements.  Three  propel¬ 
lant  types  were  used  in  this  investigation.  All  three  utilized 
a  GAP  binder  with  ammonium  perchorate  (AP)  oxidizer.  Each  had 
varied  aluminum  content,  with  the  percentages  present  being 
0%,  2%,  and  4.8%.  The  size  of  the  aluminum  powder  was  20 
microns.  The  major  ingredient  compositions  for  the  4.8% 
aluminized  propellant  were  14.7%  GAP,  45.7%  200  um  AP,  24.6% 

25  um  AP,  8.5%  TEGIN  and  4.8%  20  ym  AZ  .  The  0%  and  2% 
aluminized  propellants  kept  the  same  solids  loading. 

B.  SYSTEM  CALIBRATION 

As  previously  mentioned,  the  apparatus  was  calibrated  by 
both  Harris  [Ref.  10]  and  Kertadidjaja  [Ref.  11].  The 


apparatus  was  re-calibrated  in  this  study  as  a  result  of  the 
modifications  that  were  made.  The  standard  calibration  proce¬ 
dure  described  in  detail  by  Harris  [Ref.  10]  remained  the 
same.  The  laser  beam  passed  through  a  glass  container  holding 
a  suspension  of  polystyrene  spheres  and  water  and  was  initially 
positioned  on  the  number  one  diode.  The  linear  diode  array 
was  then  moved  vertically  down,  with  the  aid  of  a  dial  indica¬ 
tor,  such  that  the  number  one  diode  was  one  centimeter  below 
the  main  beam.  Since  the  diode  array  was  2.56  centimeters 
long  and  was  50  centimeters  from  the  focusing  lens,  this 
arrangement  allowed  for  a  range  of  measuring  angles  from  .02 
to  .07  radians.  Therefore,  the  scattered  light  striking  the 
diodes  would  not  be  affected  by  the  main  transmitted  beam  or 
the  light  diffracted  by  the  beam  stop  as  in  the  previous  two 
studies.  This  arrangement  was  used  both  for  the  exhaust  beam 
and  the  motor  cavity  beam. 

Two  polydispersion  samples  of  polystyrene  spheres,  with 
D32's  ^ . 8  microns  and  10.2  microns,  were  used  for  the  cali¬ 

bration.  These  samples  were  chosen  because  the  D^2  was  known 
and  the  expected  D^2  of  the  actual  aluminum  oxide  particles 
during  a  motor  firing  was  between  1  and  15  microns.  The 
results  are  presented  in  Figures  4.1,  4.2,  4.3  and  4.4.  These 
results  were  encouraging  in  regards  to  validating  the  modifi¬ 
cations  and  the  method  for  data  reduction.  Calibrations  were 
also  conducted  within  the  motor  cavity  with  the  laser  beam 
and  scattered  light  passing  through  the  motor  windows.  Figures 


4.1  and  4.2  are  representative  of  the  results  obtained.  This 
indicated  that  the  motor  windows  did  not  adversely  affect 
scattering  profile  and,  therefore,  bias  the  results. 

An  initial  point  of  concern  was  the  narrow  range  of  measur 
ing  angles  which  were  employed  in  the  present  apparatus. 

Since  large  particles  diffract  most  of  the  incident  light  at 
small  angles  and  very  small  particles  diffract  most  of  the 
incident  light  at  large  angles,  the  range  and  minimum  angle 
employed  could  not  have  resulted  in  sufficient  sensitivity  to 
particle  size  for  the  expected  range  of  particle  sizes.  Also 
shown  in  Figure  4.1  are  two  additional  theoretical  profiles 
for  D32  °f  10  and  8.3  microns.  It  was  apparent  from  the  sensi 
tivity  of  the  slope  change  with  D32  that  the  mean  particle 
size  could  be  determined  with  an  uncertainty  less  than  ±.5 
microns. 

The  "upsweeps"  in  the  curves  at  larger  angles  using  the 
two  angle  method  were  the  result  of  the  increase  in  slope  of 
the  scattering  profile  due  to  the  effect  of  truncation  by 
the  apparatus.  These  parts  of  the  curves  remain  only  to 
provide  a  reference  to  the  angle  ratios. 

Another  informative  point  in  apparatus  design  is  given 
by  Hirleman  [Ref.  18]  in  regards  to  placement  of  the  focusing 
lens.  Hirleman  notes  that  truncation  of  the  larger  scatter¬ 
ing  anbles,  which  is  a  function  of  the  placement  of  focusing 
lens,  will  cause  a  biasing  against  the  small  particles.  The 
expression  given  in  this  discussion  is: 


(TTd/7.66X)D 


where: 

Zc  =  critical  distance  in  front  of  the  focusing 
lens  where  significant  energy  diffracted 
by  particles  less  than  diameter  D  is 
truncated  by  the  focusing  lens 

X  =  wavelength  of  the  incident  light 

For  a  given  Zc,  all  particles  that  are  greater  than  D  in  size 
will  have  84  percent  of  their  diffracted  energy  pass  through 
the  focusing  lens.  In  the  exhaust  configuration  of  the  present 
apparatus  Zc  was  equal  to  30.5  centimeters,  d  was  5.08  centi¬ 
meters  and  X  was  .6328  microns.  This  yields  a  D  of  9  microns. 
Thus,  according  to  this  criterion,  some  of  the  scattered  light 
from  particles  below  9  microns  is  truncated  and,  therefore, 
results  in  a  bias  towards  larger  particles.  This  type  of 
biasing  was  small  in  the  present  apparatus  down  to  particle 
sizes  of  approximately  3  microns  (the  smallest  particles  in 
the  4.8  micron  polydispersion)  as  evidenced  by  the  good 
calibration  results  is  shown  in  Figure  4.3.  However,  the 
apparatus  could  be  improved  by  either  moving  the  lens 
closer  to  the  exhaust  (at  the  risk  of  apparatus  damage  from 
the  exhaust  jet)  or  by  using  a  larger  diameter  lens. 

Transmittance  effects  were  also  investigated  in  conjunc¬ 
tion  with  the  calibration.  As  previously  mentioned,  trans- 
mittances  above  90%  are  required  to  meet  single  scattering 
requirements.  Calibrations  were  conducted  using  polystyrene 


spheres  with  a  D32  of  10.2  microns.  Scattering  data  were 
recorded  for  transmittances  of  85%,  70%,  60%,  50%,  and  30%. 

The  results  are  presented  in  Figures  4.5  through  4.14,  and 
summarized  in  Table  III.  Both  the  curve  fit  and  the  two- 
angle  methods  resulted  in  very  good  agreement  with  the  actual 
mean  particle  size  with  transmittances  as  low  as  50%.  Even  a 
transmittance  of  30%  did  not  severely  affect  the  measurement. 
These  results  indicate  that  particle  sizing  methods  based  upon 
single  scattering  theory  can  effectively  be  used  in  a  multiple 
scattering  environment  with  transmittances  significantly  less 
than  90%. 

A  test  was  then  conducted  to  determine  the  transmittance 
in  the  actual  motor  exhaust  using  the  2%  aluminized  propellant. 
The  procedure  was  the  same  as  the  one  used  in  Reference  11. 

The  transmitted  beam  was  directed  onto  the  diode  array  (with 
the  appropriate  filters  to  reduce  the  intensity)  and  the  beam 
profile  was  recorded.  The  motor  was  then  fired  and  the  profile 
was  again  recorded.  The  results  are  presented  in  Figure  4.15. 
The  results  indicated  a  transmittance  greater  than  90%,  with 
some  beam  spreading.  There  was  also  no  apparent  beam  shift. 
Therefore,  transmittance  loss  was  determined  not  to  be  a 
debilitating  factor  in  the  motor  exhaust  jet. 

The  calibration  curves  have  an  apparent  limit  of  approxi¬ 
mately  .04  radians  for  which  scattered  light  can  be  recorded 
(for  example,  see  Figure  4.3).  However,  it  was  stated  that 
the  apparatus  was  designed  to  record  angles  up  to  .07  radians. 


This  is  due  to  the  water-glass  interface  having  an  effective 
index  of  refraction  of  1.39.  Thus,  the  scattering  angles  are 
corrected  for  this  high  index  of  refraction,  resulting  in  a 
corrected  maximum  angle  of  approximately  .04  radians. 

C.  NOZZLE  EXHAUST  MEASUREMENTS 

One  test  was  conducted  with  a  nonmetallized  propellant. 

The  results,  presented  in  Figure  4.16,  show  no  scattering. 
Therefore,  the  scattering  observed  by  Kertadidjaja  [Ref.  11] 
at  smaller  angles  for  a  similar  test  was  probably  the  result 
of  beam  broadening  and  diffraction  around  the  beam  stop. 

Three  tests  were  conducted  with  both  the  4.8%  and  2% 
aluminized  propellants.  The  motor  firing  data  and  a  summary 
of  the  particle  size  measurements  are  presented,  respectively, 
in  Tables  IV  and  V.  From  light  scattering  measurements  each 
propellant  showed  a  similar  result,  one  test  with  a  D.^  of 
approximately  5  microns  and  one  with  a  larger  value  of 
(7-12  microns) .  This  test-to-test  variation  probably  resulted 
from  the  very  short  burn  times  and  the  corresponding  lack  of 
a  well-defined  steady-state  condition.  These  results  are 
presented  in  Figures  4.17  through  4.28.  Typical  S.E.M.  photo¬ 
graphs  are  shown  in  Figure  4.29.  The  light  scattering  results 
were  quite  consistent  with  those  obtained  from  the  S.E.M. 
evaluation.  Although  the  sample  sizes  from  which  the  S.E.M. 
evaluations  were  based  were  relatively  small  (see  Table  VI) , 
it  still  provided  a  good  indication  as  to  the  validity  of 


the  light  scattering  measurements.  It  should  be  noted  that 

an<3  to  an  even  greater  extent  D43,  is  strongly  dominated 
by  the  larger  particles.  One  or  two  large  particles  can 
significantly  change  Future  S.E.M.  evaluations  should 

utilize  much  larger  sample  sizes. 

The  SPP  equation  for  (which  is  primarily  a  function  of 

throat  diameter)  predicted  much  smaller  values  of  D43  than 
observed  in  this  limited  set  of  data.  The  SPP  model,  however, 
is  based  on  particulate  data  from  motors  with  nozzle  diameters 
greater  than  1  inch. 

D.  NOZZLE  ENTRANCE  MEASUREMENTS 

Light  scattering  measurements  through  the  motor  at  the 
nozzle  entrance  were  unsuccessful  as  a  result  of  combustion 
light  existing  at  the  same  frequency  as  the  laser  light.  This 
was  the  case  for  both  the  helium-neon  laser  (A  =  .6328  ym) 
and  the  argon  laser  (A  =  .488  ym) . 

The  analyses  of  the  S.E.M.  data  for  the  nozzle  entrance 
indicated  a  much  larger  D^2  and  D43  as  compared  to  the  nozzle 
exit  measurements.  Therefore,  a  significant  reduction  in  size 
appears  to  have  occurred  to  the  larger  particles  across  the 
exhaust  nozzle. 
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TABLE  III 


EFFECT  OF  TRANSMITTANCE 
ON  MEASURED  PARTICLE  SIZE 
(D32  from  Manufacturers  Data  Equals  10.2) 


Transmittance 

(%) 


Measured  D. 


(microns) 


TABLE  IV 


SUMMARY  OF  MOTOR  FIRINGS 


Date  of 
Test 

Wt  Of 
Alum 
(%) 

Press 

Pc 

(psig) 

Press 
Max  P 
(psig? 

Burn 

Time 

(sec) 

Residence 

Time 

(msec) 

Oct  31,  85 

0 

350 

366 

3.8 

- 

Oct  22,  85 

4.8 

312 

328 

2.9 

10.9 

Oct  23,  85 

4.8 

328 

359 

2.8 

11.3 

Oct  18,  85 

4.8 

333 

340 

2.8 

11.3 

Nov  19,  85 

4.8 

450 

466 

.95 

14.0 

Oct  30,  85 

2.0 

296 

296 

2.8 

10.4 

Nov  26,  85 

2.0 

408 

467 

2.0 

14.4 

Nov  25,  85 

2.0 

566 

575 

1.8 

16.7 

Note:  Pressure  colume  labeled  Pc  indicates  pressure 

at  which  data  was  taken. 


TABLE  V 


SUMMARY  OF  EXPERIMENTAL  RESULTS 


Press 

Pc 

(psig) 

D32 

Light  Set. 
(psig) 

d32 

S.E.M. 

Ext./Mtr. 

D43 

S.E.M. 

Ext./Mtr. 

D43 

SPP  Model 

328 

5.6 

5.6/18.2 

7.3/24.7 

2.3 

333 

5.5 

-  /  - 

-  /  - 

2.3 

450 

12.0 

13.6/18.9 

19.7/23.9 

2.5 

296 

7.0 

12.2/14.4 

18.6/19.6 

2.3 

408 

9.3 

-  /  - 

-  /  - 

2.3 

566 

4.6 

4.9/10.8 

6.0/13.2 

2.5 

(-)  indicates  sample  was  not  taken  for  this  run. 
Diameters  are  in  microns. 


Figure  4.1  Calibration  With  10.2  Micron  Polystyrene  Spheres 


Figure  4.2  Calibration  With  10.2  Micron  Polystyrene  Spheres 


CURVE  FIT  RESULTS 
INTENSITY  vs.THETR 


3DU110A  Q3U31 


bration  With  4.8  Micron  Polystyrene  Spheres 


RNGLE  METHOD 


bration  With  4.8  Micron  Polystyrene  Spheres 


CURVE  FIT  RESULTS 


THETH  (rad) 

Figure  4.5  Curve  Fit  Method,  85%  Transmitted  Light 


THETA  (rad) 

Figure  4.7  Curve  Fit  Method,  70%  Transmitted  Light 


RNGLE  METHOD 


Two  Angle  Method,  70%  Transmitted 


THETR  (rad) 

Figure  4.9  Curve  Fit  Method,  60%  Transmitted  Light 


CURVE  FIT  RESULTS 
INTENSITY  vs.THETR 


Figure  4.11  Curve  Fit  Method,  50%  Transmitted  Light 
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urve  Fit  Method,  30%  Transmi 


DIODE  NUMBER 

igure  4.15  2%  Aluminum  Propellant,  Less  Than  90%  Attenuation 


Figure  4.16  Non-Metalized  Propellant,  22  Oct.,  Exhaust 
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CURVE  FIT  RESULTS 
INTENSITY  vs . THETfl 


THETR  (rad) 

Aluminum  Propellant,  18  Oct. ,  Exhaust 
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-ANGLE  METHOD 


8%  Aluminum  Propellant,  23  Oct. ,  Exhaust 


THETR  (rad) 

%  Aluminum  Propellant,  30  Oct.,  Exhaust 


CURVE  FIT  RESULTS 
INTENSITY  vs.THETR 


Figure  4.23  4.8%  Aluminum  Propellant,  19  Nov.,  Exhaust 
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Figure  4.24  4.8%  Aluminum  Propellant, 19  Nov.,  Exhaust 
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THETR  ( 

Figure  4.27  2%  Aluminum  Propellar 


Figure  4.28  2%  Aluminum  Propellant,  26  Nov.,  Exhaust 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  of  this  study  have  shown  that  the  measure¬ 
ment  of  D^2  in  rocket  motor  exhausts  can  be  made  using  quite 
simple  light  scattering  methods.  Apparatus  calibration 
results  and  S.E.M.  evaluations  of  the  collected  motor  exhaust 
products  support  this  conclusion. 

The  SPP  equation  model  underestimated  D43  for  this  inves¬ 
tigation.  This  indicated  that  the  model  predictions  cannot 
be  extrapolated  to  small  motors  with  throat  diameters  less 
than  one  inch.  The  light  scattering  method,  however,  can 
be  used  to  validate  the  SPP  model  for  the  larger  motors  upon 
which  it  was  based. 

In  order  to  improve  the  signal  to  noise  ratio  in  regards 
to  the  scattered  light  intensity,  it  is  recommended  that  the 
laser  beam  be  moved  closer  to  the  exit  plane  of  the  nozzle, 
where  the  particle  density  is  greater.  The  laser  beam  for 
this  study  was  7.3  centimeters  aft  of  the  nozzle  exit.  To 
resolve  the  combustion  light  problem  it  may  be  necessary  to  use 
a  different  laser,  such  as  a  CC>2  that  has  a  different  wavelength 
than  present  in  the  combustion  process.  Also,  other  binder 
systems  should  be  tested  to  determine  if  they  significantly 
alter  the  light  spectrum. 
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APPENDIX 


PROGRAM  LISTINGS 


10 

20  **#«»*>**  ROC 

30  PLOTS  RAU  DATA  »***»*»«** »**»*»***» 

40  «*«»««««  AVERAGES  »xm»*mi*x*Mx«x* 

50  !»*»**m»i»»*****»»  FILTERS  *tx***x#***x*ii*****» 

60  i *»«*»*«**»»**»»****  DETERMTNFS  HFAN  DIAMETER  «»**»** mix  *»»*»#** 

79  BT  INTERACTIVE  GRAPHICS  ****mx»**xx»xxx**» 

80  | xxxxmxxxixmxm  AND  THE  TUO-ANGt.E  METHOD  ixmwxxxmxmxix* 

99  timomimm***  Robert  Kelly  Harris  84  ixxxxxxmuxxmx** 

109  !**»»**i**m*M«»*»  Revised,  John  S  Rosa  85  ixmixmixxmm* 

119  iHHIimtHmillHmiHIIHHOmOIIHIHIHHIHilttitmillH 

120  OPTION  BASE  1 

131  FOH  /Hrdws/  Av2<1924) 

131  COM  A2<1024)  PUFFER 

149  TOM  /Canes/  THlt24),L 

150  COM  /Max/  H7,H5,Kt,Yt,X«,re,Xi 

1M  COH  /Readata/  B,P,H, 8391291, Q49l291,Zi4l2»},Yl(Bl?2)  BUFFER ,T2(8192)  B11FFE 

R 


179  COM  /Tu»/  A»1 < 19241, H, HI ,F 

189  DIM  Scans(R),X(l024) 

199  COM  /Plots/  P1»1291,P24129],P34t281,P49l201 

200  INTFGER  Graff  1:12480)  BUFFER 


398 

PRINT  *  OF  THE  COMBINATION  OF  THE  MEDIUM  AND  ITS  BOUNDARY* 

488 

PR  TNT  ** 

418 

PRINT  *  THIS  VALUE  IS  APPLIED  DIRECTLY  TO  THE  DATA* 

428 

PRINT  ** 

438 

PRINT  *  AIR=1.0 

PLEXIGLASS  6  WATER  =1.39* 

448 

PRINT  ** 

451 

PRINT  "ESTIMATE  OF  MOTOR  CAVITY  COMBUSTION  PRODUCTS  6  WINDOW  =  1.22* 

468 

PRINT  ** 

478 

PRINT  *  ESTIMATE  OF  ROCKET  EXHAUST  =  1.1  OR  1.0* 

488 

INPUT  M 

498 

PRINT  CHR4U2) 

598 

INPUT  *  ENTER  LASER  WAVELENGTH  (H?Ne=. 6328, Ar=. 488) *,L 

518 

INPUT  'ENTER  THE  DESIRED  PLOTTING  INTERVAL  f)F  DIODE  ARRAY  (2,4,6)\H 

528 

F=508 

Inn  FOCUL  LENGTH  OF  OBJECTIVE  LFNS 

521 

XUM 

522 

T1(l)=,92 

538 

Died*. 825 

Inn  DIODE  SPACING  MAY  BF  .03 

549 

FOR  1=1  TO  1023 

558 

X(I*1)=I*1 

!  CREATE  AN  ARRAY  OF  DIODF  NUMBERS 

561 

TMI+l)=ATN<(I*Diod)/(F»H)MT1(l> 

*  CREATE  AN  ARRAY  OF  THETA 

578 

NEXT  I 

588 

PRINT  CIIRK12) 

598 

CALL  Dtsplavt (Old ) 

! JUST  PRINTS  INITIAL  RFM6P <S  ON  CRT 

608 

IF  Old  1HEN  CALL  ReviPu(Avt<») 

'THESE  TWO  LINES  ALLOW  FUR 

618 

IF  Old  THEN  GOTO  Gms 

! REVIEWING  REDUCED  DATA  FILES 

628 

PEEP 

638 

INPUT  *  INPUT  TWO  FILENAMES  NOW  (NO-PART  ,PART)\Q.H,(M 

640 

INPUT*  ENTER  'V  FOR  8  SCANS 

'2'  FOR  4  SCANS*, P 

658 

7ri=*;INIFRNAL,4,r  !  USES  LEFT  DISK  DRIVE  FOR  RAW  DATA  FILES 

669 

CALL  Readata 

679 

CALL  Shift 

689 

PEEP 

6,0  Serpens  ! 

”04 

PRINT  CHR4(12) 

'CLEARS  SCREEN 

719 

P’ 8= ’VOLTAGE* 

721 

P2t="DIODF  NUMBER  * 

!  STRINGS  FDR  PLOTS 

730 

P3t--yntTAGE  YS.  DIODE* 

748 

P4*=’RAW  DATA* 

759 

M5=1Q24 

751 

M«=0 . 

760 

Xt=64 

!  SET  UP  VALUES  FOR  PLOTS 

770 

Yr=.l4 

788 

X«=4 

798 

Y«=5 

791 

X»=0 

898 

CALL  Plet(Y2(*)) 

!  DRAWS  AXES  ,FTC... 

818 

GSTORE  Graft*) 

!  STORES  GRAPHICS  DTSLAY  JUST  MADE 

828 

C'V.L  Datapl»t(B,Yl(*),H> 

•PLOT  NO-PARTICLES  DATA 

838 

PRINT  * 

(O-PART ICLES  PLOTTED* 

848 

PRINT  ** 

858 

PRINT  * 

HOW  MANY  SCANS  SEEM  TO  BE  GOOD  ?* 

864 

INPUT  J 

71 


■  ■ 
a 


WHICH  SCANS  ARE  GOOD  7. . .  ie. , . 1 ,2,4,5,7, INCLUDE  LAST  CO 
'LAST  COMA  IS  REWIRED  OR  TOO  HAVE  TO  HIT  CDNT..  TN1CE 


SSI  HIST 

MMAI 
min 

891  INPUT  ScansU) 

910  PRINT  CHRH12) 

911  PRINT  USING  */////////* 

™  *  AVERAGING  THE  SELECTED  SCANS* 

931  GRAPHICS  OFF 

941  CALL  A»erageU,Scans(*),Yl(*),A*l(*)) 

958  PRINT  USING  'g* 

961  GLOAD  Graft*)  (LOADS  GRAPHICS  ARRAY  RATHER  THAN  HASTE  TINE  RE-DRAHING 
971  GRAPHICS  ON 

981  CALL  Res»lt(A*l(«),X(*),H)  !  PLOTS  1024  ELEHENT  ARRAYS 

990  PRINT  USING  *//* 

lilS  PRINT  *  Averaqe  Intensity  N#-Par  tides*  • 

1010  ON  KEY  9  LABEL  ’AVERAGE 1*  GOTO  Screen 
1021  ON  KEY  1  LABa  ‘PLOT-PARTICLES*  GOTO  1060 

1*30  PR  IN’  *  PRESS  KEY  »  I  TO  RE-AVERAGE  OR  *  1  TO  CONTINUE* 

1040  Standby:  !  MANUAL  CALLS  THIS  INTERRUPT  DRIVEN  PROGRAMING 

1050  GOTO  Standby  !  LOOPS,  WAITING  FOR  USER  TO  DECIDE 

1060  PRINT  USING  *8* 

1070  OFF  KEY  1  (HELPS  AVOID  CONFUSION  BY  CLEARING  THAT  BOX 

1080  CALL  Datapl»t(B,Y2<*),H)  (PLOT  PARTICLES  DATA 

1090  PRINT  USING  */* 

11 0#  PRINT  *  PARTICLE  DATA  PLOTTED* 

1110  PRINT  *  FOR  A  HARD  COPY  OF  THIS  RAH  DATA  PRESS  KEY  *  6" 

1120  PRINT  *  * 

1130  PRINT  *  OR  TO  CONTINUE  PRESS  KEY  0  1* 

1140  ON  KEY  6  LABEL  ’HARD  i  RAH*  GOTO  Raw 
1151  ON  KEY  1  LABEL  'CONTINUE  *  GOTO  Select 
1160  GOTO  Standby 

1170  Ran:  CALL  Plot<Y2<») ,1 )  (  THE  ONE  (1)  IS  AN  OPTIONAL  PARAMETER 

1180  CALL  Datapl«t(B,Y2(*),H)  !  WHICH  IS  USED  TO  GET  HARD  COPIES 

1190  Select:! 

1200  OFF  KEY  6 

1210  OFF  KFY  1 

1220  PRINT  CHR4U2) 

1230  PRINT  USING  *////* 

1240  PRINT  *  HOU  HANY  SCANS  SEEN  TO  BE  GOOD  1* 

1250  INPUT  J 

12M  HWT  •  WHICH  SCANS  ARE  COQO  t...ie...l,2,4,5,7,INCLUDE 

LAST  COMMA- 
1271  INPUT  Sransd) 

1280  PRINT  CHR9112) 

1291  PRINT  USING  */////////* 

1301  GRAPHICS  OFF 

1310  PRINT  *  AVERAGING  THE  SELECTED  SCANS* 

1320  CALL  A»erage(J,Scans(»),Y2(»),A»2<*)) 

1330  PRINT  USING  *8* 

1340  P40=*AVERACED  SCANS* 


J 


i 


M 
k 


! SUBTRACTS  NO-PARTICLES  FROM  PARTICLES 


Plit  af  the  Difference  Between  Particles  and  N 

KEY  *  6  FOR  HARD  COPY’ 
KEY  I  1  TO  CONTINUE  " 


1350  CALL  PUt(A«2(»)> 

1360  GRAPHICS  ON 

1371  GSTORE  Graf(») 

1380  CALL  Res#lt(Av2U),X<»),H> 

1391  PRINT  USING  •//■ 

1401  PRINT  *  Average  Scattered  Intensity’ 

1410  ON  KEY  2  LABFl  ’AVFRAGE2*  GOTO  1080 
1420  ON  KEY  3  LABEL  ’SUBTRACT’  GOTO  1451 

M3*  PRINT  ’  HIT  KEY  I  2  TO  RE-AVERAGE  OR  KEY  »  3  TO 

CONTINUf ’ 

1440  GOTO  Standby 
1450  PRINT  USING  ’S’ 

1460  OFF  KEY  3 

1470  HAT  AvU  Av2-Avl  ! SUBTRACTS  NO-PARTICLES  FROM  PARTICLES 

1430  CLOAD  Craf(l) 

1490  CALL  Resalt(AvlO),XU),H) 

1500  PRINT  USING  ’///’ 

1510  PRINT  *  PUt  af  the  Difference  Between  Particles  and  N 

O-Partulps* 

1520  PRINT  ’  KEY  *  6  FOR  HARD  COPY’ 

1530  PRINT  ’  KEY  |  1  JO  CUNTINUE  " 

1540  ON  KEY  6  LABEL  ’  HARD  AVERAGE*  GOTO  1570 

1550  ON  KEY  I  LABEL  ’  FILTER  •  GOTO  1600 

1560  GOTO  Standbv 

1570  CALL  Pl«t(Avl (»), 1 ) 

1580  CALL  ResDlt(Avl(«),X(t),H) 

1590  PLOTTER  IS  3, ’INTERNAL* 

1600  PRINT  USING  ’//’ 

1H0  PRINT  ’  ENTER  THE  NUMBER  OF  TINES  YOU  WISH  TO  APPLY’ 

1621  ®RINT  ’  THE  DIGITAL  FILTER  FOR  SNOOTHING’ 

1630  PRINT  ’  EXAHPLF»»»»*»»»  10" 

1640  PRINT  ’  TAKES  ABOUT  1.5  MINUTES’ 

1650  INfUT  ’  ZERO  IS  ALSO  ACCEPTABLE", Fil 

1660  IF  Pi  1=0  THEN  Gatss 
1670  P40=’FIL1ERED  DATA’ 

1680  CALL  Filter(Avl(»),Fil) 

1690  PRINT  CHRIU2) 

1700  CAIL  PUt(Av|<D) 

1710  IF  Old  THEN  NAT  X=  Tl 
1720  CALL  Resvlt(AvlU)  ,XU),H) 

1730  PRINT  ’  PUt  of  the  Difference  Between  Particles  and  N 

o-ParticUs’ 

17*1  PRINT  ’  AFTER  APPLICATION  OF  A  DIGITAL  FTL 

TER’ 

1750  PRINT  ’  f(,r  ,  HARD  COPY  ’ 

1761  PRINT  ’  PREPARE  the  PLOTTER  and’ 

1770  FRINT  ’  PRFSS  KEY  0  8* 

1780  PRINT  ’  OR' 

1790  PRINT  ’  PRESS  KEY  0  9  TO  NORHALI7E* 

1800  ON  KEY  8  LAPEL  ’HARD  FILTERED’  GOTO  1830 
1811  ON  KEY  9  LABEL  *  NORHALIZE  •  GOTO  Gavss 


ENTER  THE  NUMBER  OF  TIMES  YOU  WISH  TO  APPLY’ 
THE  DIGITAL  FILTER  FOR  SMOOTHING’ 
EXAMPLE»»»»»»»»  10’ 

TAKES  ABOUT  1.5  MINUTES’ 

ZERO  IS  ALSO  ACCEPTABLE", Fil 


PUt  of  the  Difference  Between  Particles  and  N 

AFTER  APPLICATION  OF  A  DIGITAL  FTL 

for  a  HARD  COPY  ’ 
PREPARE  the  PLOTTER  and’ 
PRFSS  KEY  I  8’ 

OR' 

PRESS  KEY  0  9  TO  NORMALIZE’ 


v  v  /  v 


®GOTO  Standby 

CALI  t  The  1  is  for  hard  copy 

1849  CALL  Res«lt<A»1(*),X<»),H> 

1859  OFF  KEY  9 


1869  Gauss:  OFF  KEY  2  (KEY  1  OPTION  IB  STILl  VALID 

1974  OFF  KEY  8 

1889  ON  KEY  t  LABEL  *  FILTER  *  GOTO  1600  !  ALLOWS  ONE  TO  SMOOTH  OLD  DATA 
1871  !  WHEN  REVIEWING  REDUCED  FILES 

1871  01d=l  !  THIS  ALLOWS  FILTERING  EVEN  AFTER  ENTERING  GRAPHICS  MODE,  THF 
1892  !  DIODE  COORDINATES  AREN'T  NEEDED  ANYMORE  *SFE  LINES  NEAR  CALL  FILTER 

1999  PRINT  CHRM12) 

1910  !Av2  Array  is  Normalized  BUT  Avll*)  is  Still  Saved  For  Re-work  if  Needed 
1929  MAT  A»2=  Avt  !  TO  BEGIN  WITH  THE  ARRAY  IS  ASSUMED  TO  BE  NORMALIZED 
1938  1  AND  THE  USER  CHANGES  THIS  WITH  INTERACTIVE  GRAPHICS 

1949  i  IN  THE  SUBROUTINE  ‘Compare * 

1950  P1S=‘FHTFRED  VOLTAGE 


1969  P2$=*THETA  (rad)" 

1979  PIS=* IN  TENSITY  os, THETA" 

1989  P 44= "CURVE  FIT  RESULTS* 

1990  H5= , 079 
2999  Xt=,005 
2910  Yt=,1 
2920  X*=2 
2939  Ye=2 
'2931  X*=.  02 

2049  CALL  PIotlAollsl) 


(STRINGS  FOR  PLOTS 


(SET  UP  VALUES 
IFOR  PLOTS 


2059  CALL  Resslt(A»l(*),TK*),H)  (PLOT  OF  NORMALIZED  INTENSITY  »s.  THF TA 

2969  GSTORE  Crafts)  i  s^yfs  the  SCREEN  IMAGE 

2079  CALL  Conpare<A«1<*),H,H, Ml .Graft*)) 

2089  GRAPHICS  OFF 

2994  ON  KEY  4  1  ABEL  ‘OTHER  ARRAY*  GOTO  Cbonse 

2109  ON  KEY  5  LABEL  *  TWO-ANGLE"  GOTO  Bochele 

2129  ON  KEY  6  LABEL  ‘STORE  DATA*  GOIO  Storedata 

2130  OH  KEY  7  LABEL  •  QUIT  *  GOTO  ffoit 

2141  CALL  Mem-el  !  PRINTS  OPTIONS  ON  THE  SCREEN 

2159  GOTO  Standby 

2160  Bnchele:  1 


2170  CALL  Twmngle 

2184  CALL  Men vet 

2194  GOTO  Standby 

2244  Staredata: 1 
2210  CALL  Stare 

2220  CALL  Mem-e! 

2239  GOTO  Standby 

2249  Quit:  ! 


2259  END 


2260 

2279 

2289 

2294 

2394 

2314 


SUB  Averane(J,Scans(*),Y(*),Av!*))  I  AVERAGES  SELECTED  SCANS 

Av*  191  'INITIALIZES  THE  ARRAY  LOCAL  TO  THIS  ROUTINE 

FOR  1=1  TO  J  STEP  1 

K=(?can5(I)-l)*19244l  !  THIS  COUNTER  IS  THE  BEGINNING 

L=Scans(  11*1924  !  AND  THIS  ONE  THE  END  OF  BIOCKS 

FOR  I1=K  TO  L  !  OF  1924  IN  THE,  OVERALL  DATA 


74 


I 


rav\ 


2320 

2330 

234# 

2351 

23fc8 

237# 

2389 

2390 
240# 
2410 
242# 
243# 
2440 
245# 
246# 
24?# 
2480 
2490 
25(0 
251# 
252# 
2530 

2540 

2541 
2550 
2560 
2570 
2580 
2591 
2600 
2610 
2620 

2630 

2631 
2640 
2650 
2661 
2670 
2680 
2700 
2710 
2720 
2730 
2740 
2751 
2760 
2770 
2790 
2790 
2800 
2810 


I2=11-K+1 

A*U2)=A*(I2)*Y1I1) 
NEXT  II 
NEXT  I 

NAT  A»=  A»/IJ) 

SUBEND 

SUP  Dataol»t(B,Y(») ,H) 

CON  /Nil/  N7,H5,Xt,Yt,X*,Yii,Xi 
l  DIR  I 
LORC  4 

FOR  J=1  TO  P 
HOVE  5#0+J»50,H7-,#5 
LAPEL  J) 

LINE  TYPE  t 
K=(J-l)*1024*l 
L-T«1024 
HOVE  1,Y<K) 

FOR  I=K  TO  L  STEP  H 
X=I-K+1 
PLOT  X, Y(l) 

NEXT  I 
NEXT  I 
PE  NOP 
NB=0 
SUBEND 

SUP  ReswJt<Ia»<»),D(»>,H) 

a  ip  on 

HOVE  DU),Iav(D 

FOR  1=1  TO  1024  STEP  H 
DRAW  D(I) ,Iav(I) 

NEXT  I 
PFNUP 
SUPEND 
Snlt:  ! 

SUB  PI 9t(Y(»), OPTIONAL  Hard) 

CON  /Hi* /  H7,H5,Xt ,Yt,X«,Y«,X* 


1  TH 


IS  ALWAYS  HFTWEEN 


■B  IS  4  OR  8(THE  NUMBER  OF  SCANS) 


•EACH  SCAN 
!  SEE  NEXT  LINE 

!  HELPS  KEEP  TRACK  OF  FACH  SCAN  AS  IT  APPEARS 

! BEGINNING  OF  EACH  SCAN  AND 
•THE  END  WITHW  THE  TOTAL  BLOCK 
!  HOVE  TO  THE  FIRST  POINT 

•  THIS  GIVES  1  TO  1»24  FOR  ABSCISSA 


!THF  AVERAGE  INTENSITY  IS  PLOTTED 


•IF  THE  OPTIONAL(Hard)IS  PECEIVFD 
•  THE  FIGURE  GOES  10  THE  PLOTTER 


CON  /Plats/  Pl$l201,P2il20],P3#l2BI ,P4$(201 
N7=HAX(Y(*))»1.1 
H6=.1*H7 
G1N1T 


SELECT  NPAR 
CASF  1 

PLOTTER  IS  3, ’INTERNAL" 
CASE  2 

PLOTTER  IS  705, ’HPCL* 
END  SELFCT 
GRAPHICS  ON 
CSIZE  5,. 6 
DEG 
LOIR  0 


•  A  SCALING  VARIAPIE 
*  ANOTHER  SCALING  VARIABLE 
•G1NIT  IS  JUST  GOO#  PRACTICE  SO  YOU  KNny 
INHERE  YOU  ARE  BEGINNING 
•THIS  DETECTS  IF  THE  HARO  TOPT  TS  DESIRED 
? IT  COULD  BE  DONE  WITH  IF  THEN  LOGIC 
? BUT  IS  PRESENTED  FOR  FAMILIARIZATION 
•SINCE  IT  IS  MORE  POWERFUL  IN  CONPLEX 
•SITUATIONS 

!  TO  BE  ABLE  TO  SEE  THE  PLOT 


75 


i\* 

‘Si 


1 


• NUHBFR  THF  Y  AXIS 


LO»G  5 
MOVE  73,9*5 

LABEL  USING  'K\P4t 
HOVE  75,9(1 
LABEL  P30 
HOVE  75,20 
LABEL  P2* 

LOIR  90 
HOVE  18,60 
LAPEL  PI* 

VIEWPORT  38,125,30,85  •  SCREEN  UNITS  FOR  MARGINS 

FRAME  !  DRAWS  A  BOX 

WINDOW  Xx,H5,-H6,H7 

AXES  Xt ,H6,Xx ,0,Xn,l  ,5 

LOIR  I 

CSIZE  3,. 5 

LORG  8 

CLIP  OFF 

FOR  I=-H6  TO  H7  STEP  H7/18  ! NUHBFR  THF  Y  AXIS 

MCVE  X»,I 

LABEL  USING  **,HD.8DDD,',I 
NEXT  I 
CSTZE  3,. 6 
LOIR  90 
LORG  8 

FOP  I=Xx  TO  H5  STEP  Xt  (NUMBER  THE  X  AXIS 

HOVE  T.-H6 
LAPEL  USING  *I,X*;1 
NEXT  I 
PENUP 
SUBEND 

SUB  Co«parp(A»1  (*),H,H, HI, INTEGER  Graft*)} 

COM  /Ga»5s/  T1(*),L 
COM  /Wrdqau?/  Av2(*> 

COH  /Max/  H7,M5,Xf ,Yt,Xn,Y«,Xi 
Put:  ! 

PRINT  *  INPUT  REFERENCE  ANGLE  * 

INPUT  Tt 

FOR  N=1  TO  1024  STEP  1 
IF  T1(N))=Tt  THEN  3208 
NEXT  N 
lt=0 

FOR  I=N-5  TO  N+5 
It=It+A*l(I) 

NEXT  I 
A1=It/ll 

OFF  XEY  (GETS  RID  OF  ALL  LABaS  ON  KFYS 

ON  XEY  3  LABEL  "HENUE*  GOTO  4160 

PLOTTER  IS  3, 'INTERNAL1  !  IN  CASE  A  HARD  COPY  WAS  JUST  HADE 


(NUMBER  THE  X  AXIS 


J25U 

3260 

Hard=0 

Centerline=l 

1 

3279 

D=10. 

»  INITIAL  PATICLE  MEAN  DIAMETER  IN  MICRONS 

3280  Change:  CSIZE  4,. 6 

3290 

PRINT  USING  *//////* 

3380 

PRINT  * 

STAND  BY  FOR  CURVE* 

3310 

HOVE  8,1 

3320 

CLIP  ON 

3331 

DIM  T(l:1024) 

3340 

FOR  1=1  TO  1024  STEP  H 

3350 

T(I)=Al»EXR(-D*2»UTHI)‘2-Tt*2)»t((Hli.57«PI)/L)*2>)) 

3378 

DRAW  T1(I),T(I) 

3375 

IF  (THI)*((P1*D»H1  )/L) > >3. 1  THE*  34  4  0 

3380 

NEXT  I 

3440 

PENUP 

3450 

IF  Hard  THEN  3550 

•EXITS  ROUTINE  IF  A  HARD  COPY  NAS  JUST  MADE 

3469 

PRINT  USING  */////////* 

3470 

PRINT  * 

USE  THE  KNOB  TO  VARY  THE  PAR1ICLF 

SIZE* 

3480 

PRINT  * 

OR  HIT  KEY  »  6  FOR  HARD  COPY* 

3490 

PRINT  *  " 

3500 

PRINT  * 

KEY  0  9  ALLOWS  YVU  TO* 

3510 

PRINT  * 

SELECT  NEW  REF  ANGLE* 

3520 

PRINT  *  * 

3530 

PRINT  * 

HIT  KEY  *  3  TO  GET  OUT* 

3540 

PRINT  USING  *//////* 

3541 

OFF  KEY  0 

3550  OH  KEY  6  l APFL  'HARO  COPY*  GOTO  Hardqaess 

3560  ON  KEY  9  LAPEL  "NORMALIZE*  GOTO  Paat 


3570  ON  KNOB  .05  GOTO  Pulse 

3500  Wait:  GOTO  Wait 
3590  Poise:  PRINT  USING  *P* 

3600  Strnnlis’Nean  Dianeter  1 

3610  Strng2*=“  Hicrons* 

3620  n#ont=KNfiPX 

3630  D=DR0UND(D+Coont/15,2) 

3640  GLOAD  Graf <«> 

3650  LOIR  0 

3660  LOIC  8 

3670  HOVE  H5,.8»M7 

3680  LABEL  StrrgmVALI<D>*Strnq2l 

3690  GOTO  Change 

3700  Hardqavs5i ! 

3710  PRINT  * 

3720  PRINT  * 

3731  PRINT  * 

3740  PAUSE 

3750  Hard=l 

3760  CALL  Plot (A«2<«) , 1 ) 


! ( .05)  IS  TIME  INTERVAL  IN  WHICH 
!  PULSES  FROH  1HE  KNOfi  ARE 
!  COUNTFO  AND  THIS  NUMBER  IS 
!  USED  PY  THE  1NTEPRACTIVE 
!  GRAPHICS  TO  VARY  THE  PARTICLE 
!  SI7E  A*ND  PLOT  THE  ASSOCIATED 
!  GAUSSIAN  APPROXIMATIONS  OF 
!  SCATTERING  PROFILES 


PREPARE  the  PLOTTER* 
PRESS  CONTINUE  far" 
HARD  COPY* 

!  SO  THAT  SUPEXIT  OCCURS  AFTER  HARD  COPY 
!  1»*HARD  COPY 


CALL  Result(Av2(»),Tl(«).H) 

PRINT  USING  *?• 

LORG  8 

WOVE  H5,.8*M7 
LDIR  9 

LABEL  StrnqlttVAL$(9)lStrnq2t 
GOTO  Change 
Subexit:  ! 

SUBEND 
SUB  Menget 
PRINT  USING  *8* 

PRINT  ’  YOU  CAN  RF-NORMALIZF.* 

PR  TNT  "*»»****m*****g«*»****»»m*»m*»»tni)fg***Hxxm****!HHHi**’ 
PRINT  "  YOU  CAN  RE-AVERAGE  (New  Data  Only)* 

I  PRINT  'HHiMHHHIIIiliHHHmiHiOmHHIHimilllXHHlKX1 

PRINT  *  PRESS  KEY  *  4  TO  LOOK  AT  OTHER  DATA* 

PRINT  *  OLD  OR  NEW  * 

PRINT  ’  MOTOR  /  EXHAUST* 

PRINT  'HlilHmHIKiHMmmillKIMHtfdlXIHfXKIHIIlIHlf 

PRINT  *  PRESS  KEY  I  5  * 

PRINT  *  FOR  THE  '  TWO-ANGLE  '  METHOD* 

PRINT*  OF  PARTICIE  SI7ING* 

PRINT  'HHIIHIIIIIIHmMmiiHKMIHHXmHHXHHIHHKIiX' 

PRINT  *  TO  STORE  THE  REDUCED  DATA  PRFSS  KEY  9  6 

PRINT 

PRINT  *  TO  END  THIS  SESSION  PRESS  KEY  I  7* 

SUPETU 

SUB  Display! (Old) 


PRINT  * 

PRINT  ** 

PRINT  * 

PRINT  *  * 

ON  KEY  1  I  ABEL  * 

ON  KEY  2  LAREL  * 
Wait:  GOTO  Wait 
New:  PRINT  CHRSU2) 
PRINT  *  I 

PRINT  " 

PRINT  * 

PRINT  "* 

PRINT  * 

PRINT  ** 

PRINT  * 

PRINT  ** 

01d=8 
i  SUBEXIT 
Review:  01d=l 


TO  REDUCE  NEW  DATA  PRFSS  KEY  I  1* 

TO  REVIEW  PREVIOUSLY  REDUCED  DATA  PRESS  KEY  I  2* 

NFW  DATA*  GOTO  New 
OLD  DATA*  GOTO  Review 


PUT  THE  DISK  IN  THE  LEFT  DRIVE  IF  IT  IS  NOT  ALREADY* 


ENTER  THE  NAMES  OF  THE  FIRST  AND  SECOND  FILES,’ 


EACH  FILE  HAS  DATA  FROM  BOTH  DIODE  ARRAYS. 


FIRST  IS  NO-PARTICLES  &  NEXT  IS  PARTICLES  - Dlt,D2»* 


(THIS  VARIABLE  IS  PASSED  TO  THE  MAIN  PROGRAM  TO  INDICATE 
'THAT  THE  DATA  TO  PE  READ  IS  ONE  (t)  SCAN  A®  THAT 
(THE  AVERAGING  ROUTINES  WILL  NOT  BE  USED 


4580 

SUBEND 

4590 

SUB  Readata 

4*00 

CM  /Readaia/  B,P,H,Q3*I20I, 04*1201, Zz$l201,Yl(*)  BUFFER. Y2(»)  BUFFER 

4610 

Xyz=l 

!  FILE  PQINTFR  VARIABLE 

4620 

P=8 

!  NUHPFR  OF  SCANS  IN  MOT0P  DATA 

4630 

R 1=65536 

! NUMBER  OF  BYTES  OF  MOTOR  DATA 

4640 

IF  P=2  THEN  Xyz=4097 

(RECORD  *  HHERF  4  SCANS  DATA  BEGINS 

4650 

IF  P=2  THEN  R1 =32768 

1NUHRER  OF  BYTES  OF  4  SCAN  DATA 

4660 

IF  P=2  THEN  B=4 

•NUMBER  OF  SCANS 

4670 

PRINT  *  • 

4680 

PRINT  *  READING  DATA  FROM  FILE  ON  DISK’ 

4690 

ASSIGN  SDi<U<1  TO  Q3UZz* 

4700 

ASSIGN  8Disk2  TO  04UZzt 

(OPEN  I/O  PATHS 

4710 

ASSIGN  PB"ff1  TO  BUFFER  Yl(*> 

4720 

ASSIGN  PBuff?  TO  BUFFER  Y2t«> 

4730 

CONTROL  PDi'kl ,5;Xyz 

4740 

TRANSFER  POiskl  TO  gB«ffl;COUNT  R1 

IRFADS  NO-PARTICI  E  DATA 

4750 

WAIT  FOR  EOT  PDiskl 

i  NOT  AN  OVERLAPPING  TRANSFER 

4760 

CONTROL  RDisk2,5;Xyz 

i Xyz  SETS  DISK  FILE  POINTFR  TO 

4770 

(EITHER  MOTOR  OR  EXHAUST  DATA 

4700 

TRANSrEP  R0isk2  TO  0Bttff2;COUNT  R1 

(READS  PARTICLE  DATA 

4790 

WAIT  FOR  EOT  8Disk2 

4600 

ASSIGN  8Cisk1  TO  * 

4810 

ASSIGN  RDisk2  TO  « 

(JUST  GOOD  PRACTICE  TO  CI03E 

4820 

ASSIGN  oPtiffl  TO  » 

(I/O  PATHS 

4630 

ASSIGN  8Btiff2  TO  » 

4840 

SUBEND 

4850  SUB  Plot2(F,D(*),X) 

4860 

DIN  Title1*l251,Thetall20] 

4870 

GINIT 

4800 

IE  X=1  THEN  PLOTTER  IS  705, mL* 

4990 

DEG 

(DEGREES  for  LABEL  DIRECTION 

4900 

GRAPHICS  ON 

4910 

VIEUPORT  35,125,35,85 

4920 

Hax=10*INT<(HAX(D<*)M18)/l0) 

4930 

UINPOU  0(1, 2), DIE, 2), 0, Max 

4940 

G=(D(E,2)-D(1,2))/(E-1)*4 

(AN  X  GRID  LINE  EVERY  4th  POINT 

4950 

F=(INKE/4)-l)*4 

4960 

IF  E=8  THEN  F=4 

4970 

CLIP  D(4,2)-G,D(E,2)+2»G,0,Max 

(MAKES  GRID  UNIFORM 

4980 

GRID  G, 5,D(4,2)-G,I 

4990 

CLIP  OFF 

5000 

LORG  8 

5010 

LDIR  90 

5020 

CSIZE  4,. 5 

5030 

FOR  1=4  TO  E  STFP  4 

(PUTS  NUMBERS  ON  X  AXIS 

5040 

HOVE  D(I, 2) ,0 

•PUTS  NUMBERS  ON  X  AXIS 


m 

UK  I 

5998 

FOR  1=10  TO  Hax  STEP  10 

! NUMBER  Y  AXIS 

5190 

MOVE  B( 4,21-0,1 

5110 

LA8EL  USING  *t,K*}I 

5129 

NEXT  I 

5130 

CSIZE  6 ,.6 

5140 

Title»=*TUO-ANGLE  METHOD* 

5150 

Title1*=*For  Various  Angle  Ratios* 

!  STRINGS  FOR  PtOTS 

5160 

Size4=*SlZE  («icrons)' 

5)70 

Thta*=*THETA  (rad)* 

5180 

S#b*=*1* 

5190 

LOIR  90 

5290 

LORG  5 

5210 

6=0(4, 2)-C-(D(E, 2)-D(l, 2)1/10 

5220 

MOVE  B,M«/2 

5230 

LAPEL  Size! 

5241 

mm  o 

5250 

A=(D(E,2)*DU ,2)  >/2 

5260 

MOVE  A, -Max/4 

5270 

LABEL  Thtal 

5281 

LORG  3 

5290 

MOVE  A, -Max/4 

53  9  0 

LABEL  Sib! 

5310 

LORG  5 

5320 

HOVE  A,Maxil,15 

5330 

LAPEL  Title! 

5340 

CSIZE  4.5 

5350 

MOVE  A,Hax*l .05 

5361 

LAPEL  Titlelt 

5370 

FENUP 

5381  SUBEND 

5390  SUB  Distribntion(D(*T,Tratio,El 

5400 

LINE  TYPE  1 

5410 

CLIP  ON 

5420 

Tratio!=VAl9(Tratiii) 

!  ANGLE  RATIO 

5430 

MOVE  0(1,21,0(1, 1) 

5440 

FOR  1=1  TO  E 

! THIS  SUBROUTINE  PLOTS  THF 

5450 

DRAU  D(I,2),D(I,1) 

'PARTICLE  SIZES  TERIVE0  USING 

5460 

NEXT  I 

'VARIOUS  ANGLE  RATIOS  APPLIED 

5470 

LORG  4 

! TO  THE  DATA  OVER  A  RANGE  OF 

5430 

LINE  TYPE  1 

(ANGLES 

5*90 

CSIZE  4,. 3 

5500 

MOVE  DU-1, 2)  ,0(1-1,11+2 

5510 

LABEL  Tratiol 

5520 

PENUP 

5531  SUBEHD 

5541  SUB  Twoangie 

80 
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'•WfiHHiK*  SUBPROGRAM  ' TWOANGLE' 


J\ 


IIXHIHimtd 

'•HCHIHtMBI 

MIIHmiiOilH 

ImiHiiiKmi 


PARTICLE  SIZING  BY 
MEASURING  THE  RATIO 
OF  INTENSITY  AT 
TWO  ANGLES 


umiiHiiiHimii 

IHHHtHHKHill 


I  IlHIIIHIHHlllHHHHHHHHIlHHHUliillHHHHIHHHlHI 
I  OPTION  BASE  I 
I  COM  /Two/  A*1(1«?4),H,H1,F 

I  COM  /Gauss/  T1(1024),L  !  THE  GAUSS  IS  NOT  USED  HERE 

)  ! PUT  THE  CON  BLOCK  HAS  THETA 

I  DIM  0(2*11,2)  LAND  WAVELENGTH 

)  PRINT  CHRK12) 

I  PRINT  USING  V//////* 

I  PRINT  •  THIS  SUBPROGRAM  USES  THE  TWO-ANGLE  METHOD  DESCRIBED  BY  BUCHFIE* 

I  PRINT  *  ■ 

I  PRINT  '  TO  CALCULATE  PATICLE  SIZE  FOR  VARIOUS  ANGLE  RATIOS  AND  ANGLES.' 

I  PRINT  '  ' 

I  FRINT  •  IT  IS  HOPED  THAT  THE  CURVES  WHICH  RESULT  WILL  SHOW  • 

I  PRINT  ■  ■ 

I  PRINT  '  WHICH  SIZE  IS  THE  MOST  PROBABLE.' 

I  PRINT  *  * 

I  PRINT  'AFTER  NOTING  FROM  THE  RAW  DATA  WHICH  ANGLES  CONTAIN  THE  CFNTER  LOBF* 

I  FRINT 

PRINT  •  ENTER »»»*»»»»*»»«* THE  SMALLEST  USEABLE  ANGLE,' 

I  PR-'NT  •  THE  SMALLEST  ANGLE  RATIO, AND" 

PRINT  '  THE  STEP  SIZE  BETWEEN  ANGLE  RAIDS' 

I  pPW  •  YOU  WISH  TO  EXPLORE* 

PRINT  ■  EXAMPLEmmmmm***  .012, 1.2,. 4' 
i  INPUT  ■  ENTER  thetal, theta-ratio, theta-step', Q, A. B 

OFF  KEY 

X=0  (GRAPH  WILL  APPEAR  ON  SCREEN 

i  Pegin:  (RUNNING  CONTINUES  HERE  WHEN  A  HARD  COPY  IS  DESIRED 

(Van  Dp  H»lsT  and  Gvnprech  4  Sleepevich  eiolain  that  the  chanoe  in 
' wavelength  of  the  bean  is  accounted  for  by  dividing  by  the  indpx 
'of  refraction  of  the  nedivn.  #*§*«m»L=L/Hl  *»«**»***» 

C='L/Ml/,57/?I)*2  ! - see  page  15  of  nasa  terh  oaoer  ?15<> 

FOR  N=1  TO  1024  (to  see  this  is  a  convenient  constant 


I  IF  Q<T1(N)  THEN  5950 
NEXT  N 

FOR  Tratio=A  TO  3  STEP  B 
l  FOR  J=N  TO  1010/Tratio  STEP  10 

Thtal  =  Il  ( J) 
Thta2=Tl(J«Tratio) 

i  D"ltatheta=Thta2*2-Thtar2 

11=0 

I  12=0 


(FINDS  POSITION  OF  MINIMUM  ANLGE 

(VARIOUS  ANGLE  RATIOS 
'SETS  THE  RANGE  OF  POSSIBLE 
!  ANGLFS  (THETA  t) 

!  THETA  2 


M. 

6051 

6060 

6079 

6080 
6099 
6199 
6119 
6129 
6139 
6149 
6159 
6169 
6179 
6181 

6199  Xconp: 


™  ‘min.  \mm wtK  F“ 

NEXT  I  !0F  A  STEADIER  CALCULATION 

11=11/11 

FOR  I=INT(J»Trati«-5)  TO  INT( J»Tratio+5) 

I2=I2tA»l(I)  !  DONE  HERE  FOR  TMETA2  FUR 

NEXT  I  !  THE  GIVEN  ANGLE  RATIO 

12=12/11 

E=( J-N) /1 0+1  • THIS  IS  A  COUNTER  FOR  THF  ARRAY  CONTAINING 

•PARTICLE  SIZE  AND  THETA1  FOR  THE  GIVEN 
•ANGIE  RATIO 

IF  11(12  OR  11=9  OR  121=0  THEN  DtE,l)=0  (ALLOWANCE  FOR  IF  THE  DATA 
IF  11(12  OR  11=9  OR  I2(=9  THEN  GOTO  Xc»«p  IIS  NOT  VEIL  BEHAVED 
D(E,1)=SQR(-C/Deltatheta*LOG(I?/I1>)  ITUO  ANGIE  MEW'D 

! VALUE  OF  DIAMETER  based  tn  INTENSITY  RATIO 
•For  a  given  ANGLE  RATIO  and  ANGLE  THETA  1 

D(E,?'=Tbta1 


6290  Spar=PI»0(E,l)»N1/L  ITHIS  IS  PI»D/IAMRDA  -  THF  SIZE  PARAMFTFR 

6210  Tbar=SD3r*fhta2  .  ITHFTA  BAR  FOR  THF  LARGER  ANGLE 

6220  IF  Tbar >4  THEN  J=1919/Tratio  ITHIS  ENDS  THE  DO  LOOP  FOR  THIS 

6239  IF  E=1  AND  Than >3  THFN  6339  I  ANGLE  RATIO  SINGE  THE  GAUSSIAN 

6249  I  the  above  line  ends  all  calculation  if  IIS  NOT  VALID  UHEN  Tbar  >  3 
6259  I  the  first  element  (E=l)failed  the  test 


6269  NEXT  J 

6279  I  FIRST  TIMF  THROUGH-  D  HAS  THF  MOST  ELEMENTS  IT  HILL  HAVE 

6289  i  SO  SET  UP  THE  GRAPH  USING  D's  PRESENT  PARAMETERS 

6290  PRINT  CHPSU2) 

6309  IF  Tratii=A  THrN  CALL  Plot2(E,D(*),X) 

6319  CALL  Distrih*ti»n(DU>,Tr»ti*(E> 


6329  NEXT  Tratio 


6339  ON  KEY  2  LAPFL  "HARD  COPY*  GOTO  Hard 

6340  ON  XFY  3  LABEL  "NENUE  *  GOTO  Subeut 
6350  PRINT  USING  */////* 

6369  PRINT  *  YOU  CAN  GET  A  HARD  COPY  BY  PRFSSING  XFY  t  2' 

6379  PRINT  *  OR  EXIT  THIS  ROUTINE  BY  PRESSING  KEY  *  3’ 

6389  Standby :GOTO  Standby 

6399  Hard : X=1  !A  VARIABLE  TO  CONTROL  PLOTTER 

6409  GOTO  Begin 

6419  Subexit:  GINIT 

6429  GRAPHICS  OFF 

6431  SUBEND 

6449  SUB  Shift 

6451  COM  /Readata/  B,P,H,Q3H29J,Q4»l291lZiit291,Y1(o)  BUFFER, Y2(»>  BUFFER 
6460  DIM  Ed : 8 1 92 ) 

6479  PRINT  CHRH12) 

6489  PRINT  USING  "///////* 

6491  PRINT  *  RAN  DATA  IS  BEING  SHIFTED  TO  CORRECT  FOR  SMALL  GAPS  BETWEEN  SCANS’ 
6500  PRINT  *• 

6519  PRINT  *  BE  PATIENT,  IT'S  A  LONG  SET  OF  DO-LOOPS’ 


82 


•\V.’ 


-  -  n*.  m 


-»  ..m'. 


a 


THERE  ARE  SOME  SHALL  GAPS  BETWEEN  SCANS  AND  THIS  SUBROUTINE 
SHIFTS  THE  DATA  SO  THAT  THE  FIRST  DIODE  DATA  POINT  IS  MOVFl) 


! MOTOR  DATA  HAS  8  SCANS 


IOPPERATES  ON  NO-PARTICLE  AND  PARTICLE  SETS 


!  SHIFTS  THE  DATA  SO  THAT  THE  FIRST  DIODE  DATA  POINT  IS  HOOF D 

!  TO  THE  VERY  BEGINNING  OF  ITS  1D24  BLOCK  IN  THE  OVERALL  ARRAY 

!  .THE  FIRST  SET  IS  RIGHT  ON,  THE  NEXT  IS  ONE  OFF, THE  THIRD  IS 
•  TWO  OFF  ,,SO  FORTH.  THIS  HAY  HOT  HATTER  WITH  OUR  RESOLUTION 

!  AND  IS  PROBABLY  DUE  TO  THE  MEMORY  CARD  CYCLING  AT  THE  END  OF 

!  SOHE  SCANS.  SEE  THE  CIRCUIT  TIMING  DIAGRAH  IN  THE  THESIS. 
SELECT  B 

CASE  4  ! EXHAUST  DATA  HAS  4  SCANS 

11=3 

CASE  8  ! MOTOR  DATA  HAS  8  SCANS 

M=7 

END  SELECT 
FOR  K=1  TO  2 

IF  K=1  THEN  HAT  E=  Y1  IOPPERATES  ON  NO-PARTICLE  AND  PARTICLE  SETS 

IF  K=2  THEN  HAT  E=  Y2 
FOR  J=0  TO  H 

IF  H=7  AND  J<=3  THEN  8800  'One  4098  Block  Doesn't  Need  Shifting 
FOR  I=(J)*1024+1  TO  (J+I)»1024  ! "LOCKS  OF  1(24 

IF  H=3  THEN  L=I*J  I Array  B  Has  the  worst 

IProhlen  With  Shifting  Data 

IF  H=7  THEN  L=I*1  I  Array  D  is  Allways  off  bv  one 

IF  L)(J+1)»1024  THEN  L=(JM)»1024 
•JUST  TO  AVOID  PROGRAMMING  ERROR  AT  THE  END  OF  THE  ARRAY 
E(I)=E(L)  (THIS  SHIFTS  THE  DAIA 

(DEPENDING  ON  'l  ARRIVED  AT 
!BY  OBSERVING  RAW  DATA 

NEXT  1 
NEXT  J 

IF  lt=t  THEN  HAT  Yl=  E 
IF  K=2  THEN  HAT  Y2=  E 


6030  NEXT  K 
6R48  SUBEND 
6851  SUB  Store 


685#  SUB  Store  !  THF  LARGE  ARRAY  CONTAINING  MULTIPLE  SCANS  HAS  BEFN  RFDUCFD 

666 D  i  TO  A  MEAN  SCATTERING  PROFILE  BY  AVFRASING  AND  FILIFR1NG,  IF 

6878  !  YOU  FEFL  THAT  STORING  THIS  DATA  IS  NECESSARY  THIS  ROUTINE 

6880  i  DOES  IT.  IT  SAVES  DISK  SPACE  TO  STORE  THE  RESULTS  THEN 

6898  !  FI IHINATE  THE  RAW  DATA  IF  YOU  FEEL  CONFIDENT  THAT  THE 

6900  i  REDUCTION  IS  THE  BEST  THAT  CAN  PE.  IN  OTHER  WCRDS, 

6918  !  DO  HOT  PURGE  A  RAW  DATA  FILE  UNLESS  YOU  ARE  ABSOLIJTEI  Y  SURE 

6928  i  YOU  WON'!  WANT  TO  REDUCE  IT  AGAIN. 

6938  COH  /Hrdqaus/  Aw2(*>  !Ao2(»)  IS  THE  REDUCED  DATA 

6948  DIH  DataU  :  IP24)  BUFFER  ! BUFFER  USED  TO  TRANSFER  TO  DISK 

6958  HAT  Data=  Av2 
6960  PRINT  CHR*(12) 

6978  PRINT  USING  *//////• 

6981  PRINT  *  A  SUGGESTED  METHOD  OF  NAMING  REDUCED  DATA  FILES  IS  AS  FOLLOWS" 

6990  PRINT  *  H - MOTOR  BEAH' 

7081  PRINT  *  X - EXHAUST  BEAH’ 


Ah  i  v -a:;  ‘ 


i mm 


7118  PRINT  '  C - CALIBRATION  ,  IF  NO  'C'  THEN  AN  ACTUAL  FIRING  15  ASSUHFD 

7028  PRINT  *  G - GLASS  BEADS  1 

7031  PRINT  *  IF  NO  'C'  THEN  'C'  STANDS  FOR  'GAP'  PROPFLLANT 

7040  PRINT  1  A - ALUNINUH  OXIDE* 

7050  PRINT  •  PPP - CHAMBER  PRESSURE  FOR  RUN  OR  OTHER.. . .(NOZZLE  TYPE)’ 

7060  PRINT  *  PARTICLE  SIZE  FOR  CALIBRATION* 

7970  PRINT  *  HHDD - MONTH, DAY  OF  RUN  OR  CALIBRATION* 

7080  PRINT  ** 

7090  PRINT  *  EXAMPLE:  MCA125JN12* 

7103  PRINT  ** 

7110  PRINT  *  MOTOR  BEAM  CALIBRATION  USING  ALUMINUM  OXIDE  FROM  1  TO  25  MICRONS 
7120  PRINT  *  ON  JUNE  12" 

7130  PRINT  *  EXAMPLE:  XG550AU10* 

7140  F9INT  *" 

7150  PRINT  *  EXHAUST  BEAM  DATA  FOR  GAP  PROPELLANT  AT  550  pci  ON  AUGUST  10* 
7160  PRINT  *  PLACE  A  DISK  IN  THE  RIGHT  HAND  DRIVE* 

7170  PRINT  *  ENTER  THE  FILENAME  YOU  WISH  TO  USE  FOR  THIS  REDUCED  DATA* 

7190  INPUT  A4 

7190  CREATE  BDAT  AO.512,16 

7200  ASSIGN  RDisk  TO  Al 

7210  ASSIGN  RP»ff  TO  BUFFER  DataUO 

7228  CONTROL  8Buff,3;l ,8192,1  (BUFFER  IS  FULL 

7230  TRANSFER  PBuff  TO  §Disk;COUNT  8192 
7240  WAIT  FOR  FOT  RDisk 
7250  ASSIGN  Muff  TO  * 

7260  ASSIGN  PDjsk  TO  • 

7270  SUPEND 

7280  SUB  RevipuXAvl (») ) 

7290  DIM  OaU(t:  1024)  BUFFER 

7300  PRINT  *  EACH  REDUCED  FILE  CONTAINS  ONE  SCAN  ONLY.  YOU  MUST  KNOU  IF  IT* 
7310  PRINT  *  IS  EXHAUST  ,  MOTOR  CAVITY  ,  OR  CALIBRATION  DATA.* 

7320  PRINT  ** 

7330  PRINT  *  THE  DISK  UITH  THE  REDUCED  FILE  SHOULD  BE  IN  THE  RIGHT-HAND  DRIVE* 
73*0  PRINT  *" 

7350  PRINT  *  ENTER  THE  FILENAME  OF  THF  REDUCED  DATA* 

7360  PPINT  ’  TO  BE  REVIEWED* 

7370  PRINT  ** 

7380  INFUT  A* 

7390  ASSIGN  BPisk  TO  A$ 

7400  ASSIGN  Muff  TO  BUFFER  Data(») 

7410  CONTROL  ?Disk,5;l 

7420  CONTROL  Muff ,3;1 ,0,1  (THIS  IS  AN  EUPTY  BUFFER 

7430  TRANSFER  8Di<k  TO  BBaffjCOUNT  8192  !  1 024*8=NUMBER  OF  BYTES 

7441  WAIT  FOR  EOT  POisk 
7450  ASSIGN  Misk  TO  * 

7460  ASSIGN  Muff  TO  • 

7470  HAT  A*l=  Data 
7480  SUBEND 

7490  SUB  Filtered), Fil) 


m 

7520 

PRINT  USING  ■////////' 

7530 

PRINT  ■ 

FILTERING  THE  SCATTERING  PROFILE1 

7540 

FOR  J=1  TO  Fil 

7550 

FOR  1=1  TO  1014 

7560 

B=I+1 

•THIS  IS  A  Synetric  Howing  Averaqe 

7570 

C=I+2 

(TYPE  OF  DIGITAL  FILTER.  EACH 

7580 

D=I+3 

(DATA  POINT  IS  EQUALLY  WEIGHTED 

7590 

F=IM 

(IN  THIS  CASE  BUT  THIS  CAN  BE 

7601 

Q=I+5 

(CHANGED  IF  ONE  DETERMINES  THAT 

7610 

H=I+6 

(FEWER  POTNTS  WITH  UNEQUAL  WEIGHTS 

7620 

K=I»7 

(WOULD  BE  FASTER  OR  GIVF  BETTER 

7630 

L=I+8 

(RESULTS.  THIS  TYPE  OF  FTLTFR  WAS 

7640 

N=I+9 

(USED  SINCE  IT  INTRODUCES  NO  PHASE 

7650 

P=I+1B 

'LAG  (Angular  Resolution) . 

7660 

E(G)=(A(I)+A<8)+A(C)+A(D)*A(F)+A(G)+A(H)+A(K)+A(LMA(N)+A(P))/11 

7671 

NEXT  1 

7680 

HAT  A=  F 

7691 

NEXT  J 

7700 

SUBEND 

85 
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